The astrophysical rate of ^^0(Qf,7)^^Ne via recoil-decay 
tagging and its implications for nova nucleosynthesis 



B. Davids 



Kernfysisch Versneller Instituut, Zernikelaan 25, 9747 A A Groningen, 
- - ■ The Netherlands 

m ■ 
O ; 

O i J. Jose 

(N . 

' Dept. Fisica i Enginyeria Nuclear, Univ. Politecnica de Catalunya, and 

. Institut d'Estudis Espacials de Catalunya, Ed. Nexus-201, Gran Capita 2-4, 08034 

, Barcelona, Spain 

Abstract 

\ The "'^^0(a,7)"'^^Ne reaction is one of two known routes for breakout from 

Q ' the hot CNO cycles into the rp process. Its astrophysical rate depends 

. on the decay properties of excited states in ^^Ne lying just above the 

\ ^^O + a threshold. We have measured the a-decay branching ratios for 

ff^ ' these states using the p(^^Ne,i)^^Ne reaction at 43 MeV/u. Combining 

■ our branching ratio measurements with previous determinations of the 

>^ , radiative widths of these states, we calculate the astrophysical rate of 

^ I ^^0(a, 7)^^Ne. Using this reaction rate, we perform hydrodynamic cal- 

Q ■ culations of nova outbursts and conclude that no significant breakout 

^ . from the hot CNO cycles into the rp process occurs in classical novae 

F; ■ via 1^0(0, 7)i9Ne. 

> : 

^ , Novae are thermonuclear runaways powered by the accretion of hydrogen- 

^ \ rich material from a stellar companion onto the surface of a white dwarf in 

■ ■ ' a close binary system. Energy production and nucleosynthesis in such stellar 

sites are in general determined by the CNO cycles, with contributions from 
the NeNa and MgAl cycles PP . The possibility of breakout from the hot CNO 
cycles into the rp process has been extensively discussed in previous work, 
in particular under the high temperature and density conditions attained in 
the most massive ONe nova outbursts. Several reactions have been suggested 
as pathways for this breakout [2 , but only two are currently thought to be 
possibilities: ^^0(a,7)"'^^Ne and ^^Ne(a,p)^^Na. In astrophysical environments 
the ^^0(Q;,7)^^Ne reaction proceeds predominantly through resonances lying 
just above the ^^O + a threshold at 3.529 MeV in ^^Ne. For nova conditions 
in particular, the reaction rate is determined by the a width Tq, of the 4.033 



MeV, 3/2+ state, owing both to its close proximity to the ^^O + a threshold 
and its low centrifugal barrier to a capture. 

Direct measurements of the low energy cross section, which require high- 
intensity radioactive ^^O beams, are planned. For states in ^^Ne lying at 
excitation energies relevant to novae and accreting neutron stars, only the a- 
and 7-decay channels are open, as the proton and neutron separation energies 
are 6.4 and 11.6 MeV [Sj respectively. Hence, by populating these states and 
observing the subsequent a and 7 decays, one can deduce the branching ratio 
Bq, = r^/r. If is also known, one can then calculate and thereby the 
contribution of each state to the resonant rate of ^^0(a,7)^^Ne. Efforts of this 
kind to detect a particles from the decay of ^^Ne states populated via transfer 
reactions were made but these experiments were insufficiently sensitive 

to measure for the critical 4.033 MeV state, which was expected to be of 
order 10"^ 0. 

In an experiment at the Kernfysisch Versneller Instituut 0, we have suc- 
ceeded in obtaining branching ratio data at this level of sensitivity. Populating 
the important states via the ^^Ne(p, t)^^Ne reaction in inverse kinematics with 
a ^^Ne beam energy of 43 MeV/u, we detected either ^^Ne recoils or their ^^O a- 
decay products in coincidence with tritons in the Big-Bite Spectrometer (BBS) 
The large momentum acceptance of the BBS (Ap/p = 19%) allowed de- 
tection of either ^^Ne recoils or ^^O decay products along with tritons emitted 
backward in the center of mass system. Positioning the BBS at 0° maximized 
the yield to the 4.033 MeV, 3/2"*" state in ^^Ne. This state, whose dominant 
shell-model configuration is (s(i)^(lp)~^ [HI, was selectively populated by an 
£ = 0, two-neutron transfer from the 3/2"*" ground state of ^^Ne. Position mea- 
surements in two vertical drift chambers (VDCs) ^U] allowed reconstruction 
of the triton trajectories. Excitation energies of the ^^Ne residues were deter- 
mined from the kinetic energies and scattering angles of the triton ejectiles. 
The 7 decays of states in ^^Ne were observed as ^^Ne-triton coincidences in the 
BBS, whereas a decays were identified from ^^0-triton coincidences. 

Recoils and decay products were detected and stopped just in front of the 
VDCs by fast-plastic/slow-plastic phoswich detectors jTT] that provided energy 
loss and total energy information. A separate array of phoswich detectors was 
used to identify tritons after they passed through the VDCs. Timing relative 
to the cyclotron radio frequency signal was also employed for unambiguous 
particle identification. An Al plate prevented many of the heavy ions copiously 
produced by projectile fragmentation reactions of the ^^Ne beam in the (CH2)n 
target from reaching the VDCs. The spatial extent of the heavy-ion phoswich 
array was sufficient to guarantee 100% geometric efficiency for detection of ^^Ne 




Figure 1: ^^Ne-triton coincidences (7 decays of states in ^^Ne). The curve is 
the sum of a constant background and 8 Gaussians centered at the energies of 
known states in ^^Ne, the widths of which were determined by the experimental 
resolution of 90 keV FWHM. 



recoils and ^^O decay products for ^^Ne excitation energies < 5.5 MeV. This 
resulted largely from the forward focusing of the ^^Ne recoils, which emerged at 
angles < 0.36° for tritons with scattering angles < 4°. The low decay energies 
of the states studied limited the angular and energy spreads of the ^^O decay 
products. 

The ^^Ne excitation energy spectrum obtained from ^^Ne-triton coinci- 
dences, representing 7 decays of states in ^^Ne, is shown in Fig. 1. Its most 
prominent peak is due to the 4.033 MeV, 3/2+ state. The 4.379 MeV, 7/2+ 
state is the second most strongly populated. Contributions from other known 
states [12] are indicated. The experimental resolution of 90 keV FWHM is 
insufficient to resolve the 4.140 and 4.197 MeV states from one another; the 
4.549 and 4.600 MeV states are also unresolved. However, the astrophysically 
important 4.033 and 4.379 MeV states are well separated from the others. The 
curve shown in Fig. 1 is the sum of a constant background and 8 Gaussians 
centered at the known energies of the states, with standard deviations fixed 
by the experimental resolution of 90 keV FWHM. 

Fig. 2 shows the ^^Ne excitation energy spectrum obtained from ^^0-triton 
coincidences, corresponding to a decays of states in ^^Ne. The observed states 
are labeled by their energies; the 4.549 and 4.600 MeV states cannot be re- 
solved. A fit consisting of Gaussians plus a constant background is shown 
as well. The states below 4.549 MeV decay overwhelmingly by 7 emission, 
while the higher-lying states observed here decay preferentially by a emission, 
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Figure 2: ^^0-triton coincidences (a decays of states in ^^Ne). The curve is 
the sum of a constant background and 6 Gaussians corresponding to known 
states in ^^Ne, the widths of which were fixed by the experimental resolution. 
The ^^O + a threshold lies at 3.529 MeV. 



simply because the larger available decay energy enhances the barrier penetra- 
bility. The background represents a larger fraction of the total events in the 
^^0-triton coincidence spectrum than in the ^^Ne-triton coincidence spectrum, 
though it is still very small. This background is well reproduced by a constant, 
the value of which is determined to a precision of 7% (lo"). 

From these excitation energy spectra, we determine the a-decay branching 
ratios of six states in ^^Ne. No statistically significant evidence for a decays 
from the states at 4.033 and 4.379 MeV was observed. For these states the 
a- and 7-decay spectra were numerically integrated in 100 keV intervals cen- 
tered at the known energies of the states. These data were subjected to both 
Bayesian and classical statistical analyses to determine upper limits on the 
a-decay branching ratios at various confidence levels. The more conservative 
Bayesian analysis yields a 90% confidence level upper limit on Bq, for the 4.033 
MeV state of 4.3 x 10~^. Using independently determined radiative widths 
for these states, e.g. we compute their a widths and a very conservative 
99.73% confidence level upper limit on the astrophysical rate of ^^0(a,7)^^Ne. 
Details of this procedure can be found in Ref. [7j. 

With this ^^0(a,7)^^Ne rate, we have performed a series of hydrodynamic 
calculations of nova outbursts on 1.35 Mq ONe white dwarfs using an implicit, 
spherically symmetric, hydrodynamic code in Lagrangian formulation pp. This 
code describes the entire outburst from the onset of accretion through the 
thermonuclear runaway, including the expansion and ejection of the accreted 
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Figure 3: A sample of relevant reaction fluxes F (s~^) in logarithmic scale 
along the accreted envelope, for a 1.35 solar mass ONe nova accreting at a rate 
M = 2 X 10^^° Mq yr~^. The mass coordinate represents the mass below the 
surface relative to the total mass. From top to bottom, panels correspond to a 
time series from the early stages of the explosion up to the ejection stage, with 
a temperature at the burning shell equal to: 9 x 10''', 10^, 1.5 x 10^, 2.5 x 10^, 
and Tpeafc = 3.21 x 10^ K (top right panel), plus three panels corresponding to 
the last phases of the explosion, when the white dwarf envelope has expanded 
to a size of R^d ^ 10^, 10^°, and 10^^ cm, respectively. The base of the ejected 
shells hes at a mass coordinate of —3.05 x 10~^. 



envelope. We assume a mass accretion rate of 2 x 10^^° M0 yr^^ and 50% 
mixing between the solar-like accreted material and the outermost, ONe-rich 
shells of the underlying white dwarf. As shown in Fig. 3, confirming the 
conclusions of Refs. [T^iMj, the ^^0(a,7)^'^Ne rate is less than 0.02% of the ^^O 

decay rate throughout the whole envelope, including its hottest, innermost 
shell (i.e., the base of the envelope), which is not ejected by the explosion. 
Since the timescale of the outburst is only a few times the timescale of the 
hot CNO cycles [7], this implies that only trace amounts of CNO material 
are processed via ^^0(Q;,7)^^Ne, even in the material retained by the white 
dwarf after ejection of the accreted envelope. Moreover, as Fig. 3 illustrates, 
the contribution of ^^0(q;,7) to the production of ^^Ne is greatly exceeded by 
the contribution of ^^F(p,7). In conclusion, no significant enrichment of nova 
ejecta due to hot CNO breakout via ^^0(a,7)^^Ne is expected. 
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